A proposed testbed for detector tomography 
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Measurement is the only part of a general quantum system that has yet to be characterized experimentally in 
a complete manner Detector tomography provides a procedure for doing just this; an arbitrary measurement 
device can be fully characterized, and thus calibrated, in a systematic way without access to its components 
or its design. The result is a reconstructed POVM containing the measurement operators associated with each 
measurement outcome. We consider two detectors, a single-photon detector and a photon-number counter, and 
propose an easily realized experimental apparatus to perform detector tomography on them. We also present a 
method of visualizing the resulting measurement operators. 



I. INTRODUCTION 



II. DETECTOR TOMOGRAPHY 



> 

00 

m 
O 

o 



X 
J3 



A quantum protocol or experiment can be divided into three 
stages: preparation, processing, and measurement. Quantum 
state tomography [H and process tomography lU re- 
spectively prescribe a procedure to completely characterize 
the first two stages, and have been successfully demonstrated 
experimentally. State tomography has played a crucial role in 
identifying and visualizing novel quantum states Pro- 
cess tomography is critical for verifying the operation of quan- 
tum logic gates and characterizing decoherence processes 
Jgt]. Completing this triad, detector tomography is a pro- 
cedure for determining the POVM (positive operator-valued 
measure) set of an arbitrary quantum detector ifioll . Thus, 
without any knowledge of the inner workings of the detec- 
tor we could predict its response to any input. To the best 
of our knowledge, detector tomography has not been demon- 
strated. Here, we propose an experimental testbed capable 
of performing detector tomography on measurement devices 
acting in the Fock space of the optical mode (i.e., the photon 
number Hilbert space). Precise knowledge of the detection 
POVM set is beneficial for measurement driven quantum in- 
formation processing, such as cluster-state computing [11]. It 
is also critical for state and process tomography, which use the 
measurement as a reference. Recently, it has also been shown 
that one can perform enhanced measurements by projecting 
onto non-classical states with a detector llT2llT3ll . Without any 
need for a theoretical model of the detector, detector tomog- 
raphy can establish which states the detector projects onto, 
possibly establishing whether they are indeed non-classical. 

We will begin by introducing the general experimental and 
theoretical requirements for performing detector tomography. 
We will then describe the detectors which we aim to char- 
acterize using our proposed testbed. Following this, we will 
provide a description of the testbed we have proposed to char- 
acterize the aforementioned detectors. Finally, we outline a 
method of visualizing the POVM elements in terms of Wigner 
quasi-probability distributions. 



Detector tomography is a method of experimentally deter- 
mining the POVM associated with the detector A theoretical 
description of detector tomography was introduced by Luis 
and Sanchez-Soto in Ref. ifioll in 1999, a maximum-likelihood 
technique was applied in Ref. ifTill in 2001, ancilla-assisted 
detector tomography was considered in 2004 in Ref. lll5|]. and 
an optimal processing scheme devised in 2007 in Ref. II16I1 . 

In elementary quantum mechanics, a measurement is de- 
scribed by a Hermitian operator A which can be decomposed 
into a sum of projectors tt^ with weights Ai, where each pro- 
jector corresponds to an outcome of the measurement: 



(1) 



These projectors tt^ are often deduced by the action of the 
measuring device in the classical regime. For example, since a 
calcite crystal spatially separates the polarization components 
of a laser beam it is reasoned to act similarly at the quantum 
level on single photons (i.e., TTi = \H) {H\ and7r2 = \V) {V\, 
horizontal and vertical polarization projectors aligned with the 
axes of the calcite). This example is justified by a theoret- 
ical model, namely quantum-electrodynamics, which shows 
the classical field reduces to field operators at the quantum 
level. However, some measurements, such as Bell-state mea- 
surement ^ or single-photon detection, cannot be traced back 
to any classical device. These quantum detectors typically 
rely on an unsystematic combination of a theoretical model 
(e.g., semiconductor theory) and detector characterizations 
of limited scope (e.g., of the detector noise and efficiency). 
In contrast, detector tomography would provide a system- 
atic method, with minimal assumptions, to characterize these 
quantum detectors and, thus, predict their action on an general 
input state. 

The most general type of quantum measurement is de- 
scribed by a POVM, a set of positive-semidefinite operators 
Oa corresponding to the measurement outcomes. These gen- 
eralized measurements are uniquely quantum in the sense that 
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they have no analogue in the classical regime. For instance, 
any POVM can be implemented by coupling the measured 
system to an ancilla system and then performing projective 
measurements on the combined system fvj\; the coupling has 
the potential to create entanglement between the two systems, 
which is impossible in the classical regime. 

In some ways, detector tomography is very similar to the es- 
tablished procedure of state tomography. In state tomography, 
one begins with an ensemble of systems prepared identically 
in state p. A measurement Oa, is performed on a subset of the 
ensemble for the purpose of estimating the probability pa'. 

Pa = TliOaP)- (2) 

This is repeated for a set of measurements {Oa}, producing a 
set of estimated probabilities {p^}- Through Eq. ^ and the 
precise knowledge of the form of each measurement Oa, one 
can estimate p. Significantly, the role of Oa and p are sym- 
metric in Eq. (|2]i. This means that with a set of known input 
states {pa} one could instead estimate an unknown measure- 
ment O through Eq. (|2]i, which is the goal of detector tomog- 
raphy. 

Despite these similarities, state and detector tomography 
have some important differences. One obvious difference 
is that the former seeks to reconstruct a single operator, p, 
whereas the latter seeks to reconstruct a set of operators 
{O^}, 7 = 1, . . . ,D, where D is the number of measurement 
outcomes. Both the density matrices and the POVM elements 
must be positive-semidefinite and, hence, Hermitian. How- 
ever, density matrices have a unit trace, whereas a POVM el- 
ement does not. Instead, a POVM set {O^} must satisfy, 

D 

E^7 = t (3) 

7=1 

where 1 is the identity, ensuring that the probabilities for all 
the measurement outcomes sum to one. This has implications 
for mathematical strategies for reconstruction, such as maxi- 
mum likelihood, where Eqs. (|2|i and (O must be included as 
constraints on the reconstructions. 

Generally, in tomography one requires what is known as a 
"tomographically complete" set to reference against. In state 
tomography, this translates to a set of reference measurements 
{Oa} that span the dP — 1 dimensional Hilbert-Schmidt space 
of the density matrix (the —1 term comes from the unit trace 
constraint). In detector tomography, the reference states {pa} 
must span the Hilbert-Schmidt space of the POVM set. A 
spanning set will necessarily have at least d? elements in it. 
To determine if a particular dP sized subset of the total set 
spans the space, one first vectorizes each element in the sub- 
set, then stacks these row vectors into a matrix, and then calcu- 
lates the determinant. If the determinant is nonzero this subset 
is a spanning set and hence, the set of reference states (or ref- 
erence measurements, in the case of state tomography) is to- 
mographically complete. We shall consider an example later 
in the paper 



III. THE DETECTORS 

Before performing detector tomography, one needs to ana- 
lyze the detectors in order to identify the Hilbert space they 
function in, and then find a suitable set of input states. We 
propose a testbed for two types of related detectors: the 
avalanche photodiode (APD), and the time-multiplexing de- 
tector (TMD). 

APDs are photodiodes that use the avalanche effect (i.e., 
impact ionization leading to the exponential multiplication of 
carriers) to achieve sensitivity to single photons. These de- 
tectors have the remarkable abiUty to detect a single photon. 
Unfortunately, if more than one photon is detected this infor- 
mation is lost in the avalanche. Hence, this detector acts as a 
binary detector; the two possible detection outcomes are de- 
tection of at least one photon (CLICK) and the detection of no 
photons (NO click). For each photon impinging on the detec- 
tor there is an intrinsic efficiency of the detector 77APD that the 
photon will cause an avalanche. The probability of detecting 
at least one out of n photons is given by 

P(CLICK)„ = 1 - (1 -T^apd)", (4) 

and the probability of detecting none of the n photons is, 

P(no click)„ = (1 - 77apd)" • (5) 

Thus, there is a non-zero probability of getting a NO CLICK 
event, even with photons hitting the detector The detection 
can be approximated with the two POVM elements, 

00 

NO CLICK .• 60 = ^ (1 - ?7apd)" \n) {n\ , (6) 

n=0 

00 

CLICK : Oi = i - ^ (1 - Tjj^paT H {n\ . (7) 

These approximate POVM elements do not include effects 
such as dark counts, afterpulsing, and detector saturation. The 
contributions from dark counts and afterpulsing can cause de- 
tection events without an actual photon being present. They 
will increase the values of the Oi matrix elements. How- 
ever since dark counts are independent of the actual count 
rate they should predominantly affect the matrix elements at 
low photon numbers. Conversely, afterpulsing is dependent 
on a counting event happening before, and thus, would be 
more prominent for higher photon numbers. Detector satu- 
ration would manifest itself by a dependence of 7]apd on n. 
These effects would have to be included in a more compli- 
cated detector model if one desired to derive more realistic 
POVM elements. Previous experiments with APDs II81I indi- 
cate that these approximate POVM elements will be a reason- 
able model on which to design a testbed. 

Several schemes have been proposed to overcome the lack 
of photon number resolution with APDs lfl9l I20I I21I1 . In a 
time multiplexing detector, the pulse under investigation is 
split into several spatial and temporal bins and then detected 
with APDs ifisi I22I I23I1 . The principle of operation of our 
implementation is depicted in Figure [T] 
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FIG. 1: Schematic of a time multiplexed detector. The incoming 
pulse (I) is split at a 50:50 beam splitter. The resulting pulses (II) 
are partially delayed (III) and then split again (IV). The initial input 
pulse (I) ends up in several temporal and spatial bins (IV), which are 
then detected with APDs. 



The pulse under investigation (I) is split at a beamsplitter 
and one of the two spatial modes (II) is delayed with respect 
to the other one (III). The two spatial modes are then re- 
combined at a beamsplitter, resulting in two temporal modes 
in each of the spatial outputs of the beamsplitter (IV). Addi- 
tional delay stages can be added to further increase the number 
of temporal modes. We propose to perform detector tomog- 
raphy on a TMD with two delay stages, resulting in a total 
of eight output bins. The advantage of splitting the incom- 
ing pulse into many bins is that if the pulse contains several 
photons these may be distributed into different bins and, thus, 
counted as individual photons by APDs. The measurement 
result of the TMD is the number of CLICK events summed 



over all the bins. Thus, for our detector, with two stages, there 
are nine outcomes, to 8 CLICKS. Unfortunately, k CLICKS 
does not imply there were only k photons at the TMD input. 
Photons may end up in the same bin and, thus, result in fewer 
clicks than the number of photons in the input pulse. Pho- 
tons also suffer from losses in the TMD and the detection ef- 
ficiency of the APDs. While this introduces an uncertainty in 
the measurement outcome on a shot-to-shot basis, for an en- 
semble measurement of the same quantum state the resulting 
click-statistics grasped in the vector p can be related to the 
photon number statistics as a vector a of the input state by 



P 



C Lcr, 



(8) 



where L represents the losses in the system and C is called the 
convolution matrix, which accounts for several photons end- 
ing up in the same bin after splitting. The loss matrix can be 
calculated by combining all losses in the system and modeling 
them as a single bea m sp litter with reflectivity ?/ioss at the front 
of the fiber network ll24ll25ll as. 



I otherwise. 



(9) 



with n^ni £ 0, . . . ,N. This simply describes the probability 
of nf out of n photons being transmitted. For the convolution 
matrix one has to take into account all possible routes a pho- 
ton could take through the fiber network ||26[l and the matrix 
elements can be calculated as follows; 



'0, 

Y.bPblPb2 ■■■Pb^-, 



n 



(»-E- 



Ci- 

■Pbt 



k > n, 
k ~ n, 
k < n. 



(10) 



In reality, the beamsplitters used in the TMD are never 
exactly 50% reflective due to variation in their manufacture. 
For a iV-bin TMD, pi, . . . , pisi are the probabilities of a sin- 
gle photon exiting from the fiber network in a particular bin. 
While these probabilities would all be equal in the case of 
50/50 beamsplitters, we set them by calibration measurements 
using an intense laser input pulse. The first two lines of Eq. 
( [Tol l correspond to the straightforward cases of more photons 
being detected than entering the detector (for which the prob- 
ability is zero, if dark counts and afterpulsing are negligi- 
ble) and detecting all incoming photons respectively. In the 
third case, k < n, all possible combinations of distributing 
these n photons into the bins have to be taken into account. 
These different combinations are represented by the fc-tuples 
b = {bi,b2, ■■■ ,bk) with bi & 1, . . . , and 5i ^ &2 ^ 
■ ■ ■ ^ bk. Some of the bins can be occupied by more than one 
photon. All possible distributions of the photons across the 



bins are described by the possible partitions of n into k parts. 
These can be represented by fc-tuples c = (ci, C2, • • • ,Ck) 
with Cj e {1, . . . , n}, and ci ;^ C2 ^ • • • ^ Cfc, J^k = n, 
and the definition cq = 0. For each bin, one has to consider 
all possible ways in which the cj photons ending up in this 
particular bin can be choosen. With the photons remaining 
to be distributed given by rt — J^tZo binomial coeffi- 

cient in Eq. ( fTOl ) accounts for this. However, if some of the 
elements of the tuple c are equal, the binomial coefficient will 
overestimate the number of distributions. Bins occupied with 
the same number of photons are not distinguishable, but are 
counted separately. This must be corrected for by dividing 
by the number of permutations these bins can form, i.e., the 
factorial of the number of bins with equal number of pho- 
tons. We denote this by id(c)!. For example for the tuple 
c= (4; 4; 2; 2; 2; 1), id(c)! =2! -3!. 

For the eight bin TMD that we plan to use in the pro- 
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posed detector tomography scheme, the convolution matrix 
was calculated using classical measurements of the probabili- 
ties pi, . . . ,pn and reads: 



C = 



/I 
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(11) 

The matrices developed above actually contain the POVM 
elements describing the TMD. We expect these POVM ele- 
ments to be diagonal in the photon number basis since the 
TMD has no phase reference and, thus, no sensitivity to off- 
diagonal coherences. In Eq. ([8]l, the matrix C • L represents 
the reaction of the detector to different numbers of incoming 
photons. When the matrix CL acts on a photon number distri- 
bution cr, the probability of a particular measurement outcome 
is found by taking an inner product of a with the correspond- 
ing row of the matrix. This action is similar to the action of 



the trace in Eq. (O, as we will demonstrate explicitly. This 
allows us to identify the rows of C • L as the diagonals of the 
TMD's POVM: 



0,= 



■[CL] 





















[C-L] 



(12) 



where j is the j-CLlCKS outcome, and [C • L]^ ^ is the value 
in the jth row and ith column of C • L. The density matrix of 
the incoming quantum state is related to a via. 
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(13) 



expressed in the number basis, where the ith element of a is 
ai . We leave the off-diagonal elements of p as undetermined, 
since they do not contribute after the trace. Using this form 
we evaluate the trace in Eq. (|2]i: 



Tr(6,/5) = Tr 
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(14) 



(15) 



(16) 



i=0 



This last line is equivalent to Eq. ([8]l, as we expected. 



As we have shown, one can derive a theoretical POVM set 
to describe a time-multiplexed detector However, as with the 
APD, this is a simplified description of the POVM and does 
not include afterpulsing (which can result in a click from one 
time bin triggering a click in the next), dark counts in the 
APDs at the outputs, as well as imperfections in the count- 
ing electronics (which can be considerably complicated). 

In summary, the APD and the TMD perform their measure- 
ments in the Fock space of the optical mode. Furthermore, 



both detectors have POVMs that are expected to be diagonal 
in the photon-number basis. Since this confines their action to 
a subspace of the Fock space, it should allow us to simplify 
the tomography. For the APDs all photon numbers apart from 
zero photons are combined into a single detection event and, 
to a first approximation, the only detector parameter is the ef- 
ficiency. This makes the APD a rather simple detector and a 
theoretical model of their POVM elements was easily calcu- 
lated. Thus, APDs would be the first detector with which we 
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propose to demonstrate a proof of principle detector tomogra- 
phy. The TMD, on the other hand, is non-trivial to describe, 
with both the loss and the convolution matrix having com- 
plicated forms. Incorporating realistic APDs into this model 
would be challenging, if not impossible. Consequently, the 
TMD is the second detector that we propose to perform detec- 
tor tomography on. It is a suitable candidate for a serious test 
of the detector tomography procedure. We propose to com- 
pare the reconstructed POVM set to the theoretical ones de- 
rived above to check how accurate our detector models really 
are. 



IV. THE PROPOSED TESTBED 

To perform tomography of these two detectors, we need 
probe states in the Fock space. We must also identify a set of 
states that span the space of the detector operators, i.e., a to- 
mographically complete set. Since neither detector possesses 
a phase-reference (such as a local oscillator), we make the 
safe assumption for our proposed testbed that the off-diagonal 
components of {O^} for each detector are zero in the photon- 
number basis. It follows that each O-y contains d free param- 
eters, one for each diagonal element, where d is the the di- 
mension of the Hilbert space. The tomographic ally complete 
set will then contain at least d input states. Unfortunately, the 
photon number representation of the Hilbert space is infinite, 
and so we must truncate it at some point for our mathematical 
reconstruction. This truncation sets the d we shall use in the 
tomography testbed we propose. A good point to set this trun- 
cation at is in the region in which the detector behavior has 
saturated to a constant, e.g. We expect the TMD's response 
to 100 photons and 101 photons will be nearly identical (out- 
come 8-CLlCKS with w 100% probability), making d = 100 a 
good place to truncate. Since we do not require sensitivity to 
inter-photon coherences a naturally tomographically complete 
set of input states {pa} for our proposed testbed is the Fock 
states. Each input Fock state would allow us to determine the 
corresponding diagonal element in each of the TMD's mea- 
surement operators {O^}. However, Fock states are difficult 
to produce with high fidelity. Fortunately, another candidate, 
the coherent state, is straightforward to produce with high- 
fidelity. 

Apart from being the quantum state that most resembles a 
classical state, coherent states have a unique property under 
attenuation: in the photon-number basis, a coherent state \a) 
can be written. 



i;) = ^ exp 

n=0 



2 



n), 



(17) 



which results in a Poissonian photon-number distribution. 
With attenuation 1 — r/, the coherent state transforms according 
to, 



U,,\a) = \r]a) = ^ cxp 



n=0 



-|n). (18) 



This is another coherent state with a lower mean number of 
photons. This property is important for the creation of known 
states in the low photon number regime, a necessity for to- 
mography on single photon detectors. Specifically, the high- 
fidelity coherent states produced by a laser could be easily 
attenuated to power levels suitable for the tomography while 
still retaining their form. Attenuation also reduces the impact 
of technical noise in the laser. The fractional uncertainty in 
the photon number of a coherent state, An/(n), is - once it is 
attenuated to a low photon level - much higher than the pulse- 
to-pulse jitter of the pulse energy of our laser source (typi- 
cal values around 2%). This suggests that the use of coherent 
states would also be reasonably robust against technical noise. 

It is well known, that a continuous set of coherent probe 
state vectors \a) form a basis for the Fock space. 



1^) 



da\a){a\ip). 



(19) 



However it is not as clear that they form a basis for the Hilbert- 
Schmidt space. The proof that they do lies in the existence of 
the P-function P{a): 



P{a)\a){a\d'^a. 



(20) 



However, we expect to use a discrete subset of coherent states 
\a)i rather than the full continuous basis for our proposed de- 
tector tomography testbed. We now attempt to find if this re- 
duced set is tomographically complete. We vectorize d co- 
herent states density matrices \ai) {ai \ ,i ~ 1, . . . ,d, keeping 
only the diagonals (since our detector POVM set is contained 
within this subspace). Stacking these row vectors into a ma- 
trix, we find. 



'\a,)\' i(ii«i)i' ... 



mad) 



\{d\a,W 



(21) 



The next step, finding the determinant, is difficult to do in 
general. Instead we consider a specific set of d ~ 10 refer- 
ence coherent states, in evenly spaced steps from \ai p = 1 to 



|aio| 



10. The determinant 



det 



0.37 
0.14 
0.05 
0.02 
0.01 
0.00 
0.00 
0.00 
0.00 
0.00 



0.37 
0.27 
0.15 
0.07 
0.03 
0.01 
0.01 
0.00 
0.00 
0.00 



0.18 
0.27 
0.22 
0.15 
0.08 
0.04 
0.02 
0.01 
0.01 
0.00 



0.06 
0.18 
0.22 
0.20 
0.14 
0.09 
0.05 
0.03 
0.02 
0.01 



0.02 
0.09 
0.17 
0.20 
0.18 
0.13 
0.09 
0.06 
0.03 
0.02 



0.00 
0.04 
0.10 
0.16 
0.18 
0.16 
0.13 
0.09 
0.06 
0.04 



0.00 
0.01 
0.05 
0.10 
0.15 
0.16 
0.15 
0.12 
0.09 
0.06 



0.00 
0.00 
0.02 
0.06 
0.10 
0.14 
0.15 
0.14 
0.12 
0.09 



0.00 
0.00 
0.01 
0.03 
0.07 
0.10 
0.13 
0.14 
0.13 
0.11 



0.00" 
0.00 
0.00 
0.01 
0.04 
0.07 
0.10 
0.12 
0.13 
0.13 



^0 



(22) 

implying that this set is tomographically complete, for a space 
truncated at photon number basis state n ~ d. Consequently, 
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it appears that a finite set of coherent states are a good set of 
reference states for our proposed detector tomography testbed. 

An essential requirement for a tomography system is pre- 
cise and accurate knowledge of the probe states. In the con- 
text of our proposed testbed, this translates to knowing a ac- 
curately. Keeping in mind that will be less than 100, a 
direct measurement is out of reach of commercial power me- 
ters. For a coherent state the mean number of photons (n) is 
connected to a via. 



(23) 



For pulsed light of the wavelength A and the repetition fre- 
quency / the time averaged power P can be calculated by. 



P 



{n)hcf 
A ' 



(24) 



which then enables a calculation of jap from the measured 
power: 



PA 



(25) 



While A and / can generally be determined accurately, P can- 
not. Systematic error in the power measurement would re- 
sult in a global scaling of the detector response: the recon- 
structed "efficiency" of the detector would scale by the er- 
ror, while the form of the POVM element would remain un- 
changed. We expect this |ap systematic error to be less than 
5%. For lap = 1, A = 800 nm and / = 100 kHz (experi- 
mental parameters, we plan to run the proposed tomography 
on), P = 25 atto Watts, a strikingly low average power. How- 
ever, since coherent states are invariant under attenuation, a 
highly transmissive beamsplitter could be used to pick off a 
large portion of the incoming beam. If the ratio of power in 
the transmitted and reflected arms is well calibrated, the high 
power arm can be used to monitor the power in the low power 
arm. This can leverage the power into the micro Watt range, 
accessible to power meters with calibrations that can be traced 
to a National Standards Institute. In our proposed testbed, this 
beamsplitter would be placed after a variable attenuator, used 
to set a, as depicted in Figure[2l For phase-sensitive detectors, 
a tomographically complete set of input states could be gener- 
ated simply by adding a variable path length for the coherent 
state in the testbed, to phase shift a, i.e., a i-^ ae"^. 

We simulated the tomography of a TMD with our proposed 
testbed by using the model given in Eq. ([8]). We used the con- 
volution matrix for our TMD displayed in Eq. ( fTTT i and simu- 
lated a detector tomography using 400 different values of jap 
for the reference states, ranging from jap = to jap ~ 40. 
In Figure[3] we plot the probabilities of both the 1 -CLICK and 
5-CLlCK outcomes against the |ap of the probe state. The 
curves for the 1 -CLICK and the 5-CLlCK outcomes have very 
different shapes with small overlap: while the 1 -CLICK ex- 
hibits a well defined peak for low jap, the 5-CLlCK peak is 
much broader and has a maximum at higher |ap. This shows 
that different measurement outcomes are associated with eas- 
ily distinguishable responses to a spanning set of input a, 
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FIG. 2: Schematic of our detector testbed. A probe beam undergoes 
variable attenuation realized with a half waveplate (HWP) and a po- 
larizing beamsplitter (PBS). A large portion of the beam is split of 
at a beamsplitter (BS) and detected by a power meter (PM) to moni- 
tor the variable attenuation, while the smaller portion of the beam is 
further attenuated with neutral density filters (ND) and then coupled 
into the detector under investigation. 




FIG. 3: Simulation of a detector tomography for the TMD. 400 dif- 
ferent values of \af between \a\^ = and \a\^ = 40 were simu- 
lated. A comparison between the 1-CLICK and the 5CLICK element 
shows a small overlap and smaller peak for the 1-CLICK element, 
while the 5-CLICK element is broader and peaked at higher photon 
number. 



suggesting that the mathematical inversion from the estimated 
probabilities (e.g., with Eq. Q) is practical. 

An analogous tomography device exists in the area of Ho- 
modyne State Tomography: the eight-port homodyne sys- 
tem. This device projects the unknown state onto a coherent 
state, much the same as our proposed testbed. It has been 
pointed out that reconstruction of the density matrix from the 
measured data in an eight-port homodyne system is an ill- 
conditioned problem | |27|] . Noise in the measured data causes 
the inversion of Eq. (O to have large errors, indicating that 
while coherent states form a tomographically complete ba- 
sis, they are still deficient for tomography. Considering noise 
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from counting statistics, in detector tomography the data set 
size is limited only by the repetition rate of the laser, which 
can be as high as 80 MHz. Thus, we can quickly accumu- 
late enough data such that counting errors for each outcome 
are insignificant. However, sources of noise, other than from 
counting, can hamper inversion. To counter the effect of these 
a common strategy is to introduce some amount of regulariza- 
tion to the inversion I28i] . Often, this regularization is implicit 



in the reconstruction, as in the Radon Transformation or Pat- 
tern Functions used in two-port homodyne tomography, which 
introduce smoothing of the data ll29ll . Indeed, filtering of data 
is a common technique established in the origins of tomogra- 
phy in medical computer imaging. These techniques have yet 
to be used in detector tomography, but we expect their appli- 
cation will be fruitful. 



V. A PHASE-SPACE REPRESENTATION OF THE POVM 

A common complaint about quantum tomography is that the end result is difficult to interpret physically IH [3^ [HI] . For 
example, in process tomography the reconstructed map for the process has x (P elements for a d dimensional Hilbert space and 
thus quickly becomes too large to understand upon inspection. A reconstructed POVM set can similarly contain a large number 
of independent parameters, (D — 1) x d^ if there are D measurement outcomes. However, the measurement operators {O^} 
benefit from their mathematical similarities to density states, allowing us to apply many of the same representation methods. 
Our detectors operate in the Fock space, which can be represented in the x and p basis, where x and p are the normalized electric 
field quadratures. This points to a particularly appealing representation commonly used to represent the state of a field mode, 
the Wigner function. 



a) Wo(x,p) 

o.5^ 



c) Wi(x,p) 




e) W5(x,p) 

0.1 



-2 -2 
-4 -4 



FIG. 4: Wigner function representations of the a) 0-CLICK, c) 1-CLICK, and e) 5-CLICKS measurement outcomes of our eight-bin time- 
multiplexed detector. A cross-section of each Wigner function is shown in b), d), and f), respectively. The Wigner functions are derived from 
a theoretical model of the POVMs of our eight-bin time-multiplexed (TMD) detector 



The Wigner function Wj is calculated in the standard way Since the measurement operators {O-,} do not have unit trace, 
from the POVM element O-y H : this detector Wigner function is not normahzed. 



dx / dpW.y{x,p) < 1, (27) 

-OO J —OO 



1 f 

W^(x,p)^— dy{x-y\d.y\x + y)e'^'Py^''. (26) .. • , u ,j u • a u u-,- a- 

' ttH J_ao ai n a I and Its marginals should not be interpreted as probability dis- 
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tributions. Nonetheless, it retains some appealing features. 
For example, the probability of measurement outcome 7 is, 

/oo /"OO 
dx / dpW{x,p)Wj{x,p), 
-00 J —oo 

(28) 

where W is the standard Wigner function of the input state p. 
Thus one can visualize the measurement as the overlap of the 
two Wigner functions. 

We plot the Wigner functions of the theoretical measure- 
ment operators of the TMD in Figure m The TMD has no 
phase sensitivity and so the Wigner function for each mea- 
surement operator is rotationally symmetric around a verti- 
cal axis through the origin. Thus, on the right of Figure |4] 
we also present a cross-section of each Wigner function to 
show its form more clearly. The Wigner functions are found 
from the TMD POVM that was derived earlier, assuming a 
loss of 48%. We show three of the nine measurement opera- 
tors, namely 0-CLlCK, 1 -CLICK and 5-CLlCK (in the figure, 
these are Wo[x,p), Wi{x,p) and, W^[x,p), respectively). 
Remarkably, the 0-CLlCK and 1 -CLICK Wigner functions look 
very similar to their state counterparts, the vacuum state and 
single-photon state. This is despite the fact that there are con- 
tributions from all incoming photon numbers to both (i.e., 
four incoming photons might all be lost and thus result in 
a 0-CLlCK event, or two incoming photons might end up in 
the same time bin and thus result in 1 -CLICK). The 0-CLlCK 
Wigner function has a Gaussian profile similar to the vacuum 
Wigner function, whereas the 1 -CLICK event goes negative at 
the origin just as the single-photon state does. 

The Wigner functions for higher click numbers quickly di- 
verge from their photonic equivalents though. The 5-CLlCK 
Wigner function is significantly different from the five photon 
Wigner function; the former has nine radial nodes, whereas 
the latter has five (the number of nodes equals the photon 
number for Fock states). Since there are only eight bins, there 
is a sizeable probability that six or seven incoming photons 
entered only five bins in total. Consequently, there will be 
significant contributions to the 5-CLlCK Wigner function from 
these higher photon numbers, distorting it from the ideal five 
photon Wigner function. In contrast, the probability that two 
or three incoming photons entered only one bin is relatively 
small, which explains why 1 -CLICK Wigner function contains 
much less distortion. 



VI. CONCLUSION 

We have detailed a proposal for performing detector to- 
mography on two detectors operating in the Fock space. The 
proposed testbed generates a sequence of calibrated weak co- 
herent states against which to reference a detector. We have 
shown that these coherent states form a tomographically com- 
plete set for the two detectors, and thus, a tomographic recon- 
struction of all the POVM elements should be possible. In the 
near future, we plan to perform detector tomography on these 
two detectors. However, we envision that a reconstruction of 
the POVM elements in the high-dimensional photon number 
space will present substantial computational challenges. The 
reconstruction will also have to deal with noise and uncer- 
tainty in the input states. This is usually not considered in 
state tomography, where the measurements are considered to 
be close to ideal. 

With the demonstration of detector tomography, experi- 
mentalists will have general procedures for characterizing all 
the parts of a general quantum device: the input states, the 
quantum circuit, and, now, the measurement. We expect this 
detector tomography to be particularly useful for devices that 
depend on generalized measurements for optimal function- 
ing, such as state or process discriminators (i.e., unambiguous 
state discrimination ifssl [34ll ). Detector tomography should 
also be useful for devices where measurement drives logic 
such as in cluster-state computing or linear optics quantum 
computing. 
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